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Abstract 

We study the production of a W boson in association with n hard QCD jets (for n > 2), 
with a particular emphasis on results relevant for the Large Hadron Collider (7 TeV and 
8 TeV). We present predictions for this process from High Energy Jets, a framework for 
all-order resummation of the dominant contributions from wide-angle QCD emissions. We 
first compare predictions against recent ATLAS data and then shift focus to obscrvables and 
regions of phase space where effects beyond NLO are expected to be large. 
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1 Introduction 

There is already a wealth of analyses of data from the Large Hadron Collider (LHC), many of 
which are beginning to stress-test the perturbative descriptions of proton-proton collisions at 
LHC energies. It is of course essential to demonstrate to what level the Standard Model (SM) 
processes are understood. This is often the first step towards a discovery (or limit) on new 
physics. 

A particularly interesting process is W plus jets where, in the leptonic decay mode, the 
neutrino leads to events with missing transverse energy. The combination of multiple jets and 
missing energy is a very common signal in models beyond the Standard Model, arising from 
a production mechanism with QCD charged particles followed by decay chains ending with a 
(semi-)stable, weakly interacting particle. 

There have been a number of data analyses of this channel already by both the ATLAS [TJ[2] 
and CMS [3] experiments. The production of a W boson in association with jets is also a clean 
environment for QCD studies, and as such could be an important testing ground for experimental 
methods used to study methods like e.g. jet vetos. What is learned in W plus jets could be 
applied in searches for the Higgs boson and other new physics. 

There has also been a great deal of recent theoretical progress in fixed-order calculations of 
this channel, where the state-of-the-art is now W plus four jets in the leading colour approxi- 
mation [4], or W plus three jets in the full colour dressing [S]-[7]. The NLO calculations lead 
to very good agreement for total cross sections and, for the right scale choices, also for several 
differential distributions. 

Much work has been reported on estimating effects beyond NLO. The NLO calculations 
for VF-production in association with up to two [8] and three jets [9] have been systematically 



merged with a parton shower. Furthermore, LoopSim 10 has also been applied to approximate 
higher fixed-order contributions to W production in regions where the X-factor is large. 

In this study we present the application of the High Energy Jets (HEJ) formalism developed in 
Refs. [ITj[l2] to the process of IF-production (and leptonic decay) in association with at least two 
hard jets. The general formalism has previously been applied to pure multi-jet production |13|, 



with the results from the accompanying flexible Monte Carlo implementation used in several 



analyses of LHC data 14-16 . The all-order results of HEJ are complementary to those obtained 
from a parton shower approach. The simplifications made by HEJ to the perturbative series to 
allow all-order results to be obtained become exact in the limit of large invariant mass between all 
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particles. This corresponds to the limit where also the (B)FKL amplitudes |17f[19| are exact. The 
production of VF+dijets was previously studied within the BFKL formalism |20| . Compared to 
this, the present formalism introduces several major improvements: it relaxes several kinematic 
assumptions made in the earlier study, improves the accuracy of the resummation, and introduces 
matching to fixed order results. 

In pure HEJ, there are no collinear singularities, no shower and no hadronisation: the output 
is a pure partonic calculation. However, this can all be consistently included by merging the 



results with a shower Monte Carlo by a procedure which avoids double-counting 21 . Such 



double-counting would otherwise arise in particular in the soft regions, which are treated in 
both HEJ and a parton shower. 

This paper is structured as follows. In Section [2] we study the various production channels 
for W plus jets at the LHC and properties of the jet radiation. In Section [3] we summarise the 
construction of the HEJ amplitudes and their implementation in a fully flexible Monte Carlo. 
The resulting program is available at HTTP:/ /cern. oh/he j. In Section [4] we compare predic- 
tions for kinematic distributions to LHC data, before embarking in Section [5] on a discussion of 
regions of phase space and observables where effects beyond NLO are expected to be large. We 
end in Section [6] with a brief discussion. 



2 The Anatomy of W Plus Jets at the LHC 

In this section, we focus on the production of a PF + /~-boson (followed by a leptonic decay) in 
association with at least two jets. We will from now denote both W + and W~ as W, and both 
electron and positron as e. Furthermore, we will apply the following acceptance cuts: 

pxj > 30 GeV, \rjj\ < 4.4, P±e > 20 GeV, |r/ e | < 2.5, Xr > 25 GeV, 
M ±w = s/2 p ±e P±u (1 - cos(c/> e - <j> v )) > 40 GeV, 

where r\ is the pseudorapidity and <j) is the azimuthal angle of the respective particle momentum. 

The kinematic requirements of the two jets and the charged lepton alone ensure that the 
parton density functions are probed only for a light-cone momentum fraction x larger than 
3.4- 10 -4 at 7 TeV (1.7- 10 -4 at 14 TeV). The neutrino momentum will contribute further to the 
momentum fraction, and thus the dynamics is well within the regime where standard, collinearly 
factorised pdfs and hard scattering matrix elements describe the cross sections accurately. We 
will therefore concentrate on perturbative higher order corrections within this framework, and 
not discuss small- x issues like /cj_-factorisation or unintegrated pdfs. We will use the standard 
MSTW2008NLO |22| pdf set throughout. The accompanying program is interfaced to LHAPDF, 
so any publicly available pdf set can be used. 



In earlier studies 13,21 23 we established a strong connection between the average number 
of hard jets in events, and the rapidity difference between the most forward/backward hard jet 
in pure jet production. This is a result not just of the opening of phase space for the emission 
of additional jets, but also the possibility of a colour octet exchange between particles ordered 



in rapidity. Such a correlation was recently confirmed by data 14 . 

It is obviously relevant to discuss the extent to which the behaviour observed in pure jets 
is relevant for the production of VF+dijets. The explanation from BFKL 17-19 24] f or the 



correlation between the average number of jets and the rapidity span of the events 23 relies 
on the higher order corrections to the dijet processes which allow for a colour octet exchange 
between all rapidity-ordered partons. 
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Subprocess 


(J 2 j [pb] 


qq' — > Wqq' 


11.6 


qg -> Wqg 


67.7 


qq — > Wqq 


9.9 


qq -> Wgg, Wq'q' 


5.8 


99 -> Wqq 


4.3 



Table 1: Leading order cross sections (in pb) for the production of a W boson in association with 
two jets for the LHC at y/s = 7 TeV. The subprocesses above the horizontal line are those which 
allow a colour octet exchange between the two initial state (and the two final state) partons, 
while those below do not. 

In Table[l]we have listed the sub-process contributions to the cross section for Wj ' j-production 
at leading order within the cuts of Eq. ([!]). The processes above the horizontal line are those 
which allow colour octet exchange between the two initial (and the two final state) partons. 
Clearly, these dominate the cross section. Such colour exchanges will dominate the higher order 
corrections in the limit of large rapidity separation between all partons. However, this require- 
ment also means that quarks are produced only as the partons which are extremal in rapidity 
in the scattering. One could worry that the cut on the centrality of the charged lepton from 
the W decay would force the quark-line emitting the W to also be central, thereby suppressing 
the phase space for a large rapidity span between the jets. In figure[l]we plot l/ada/d(yf — yi) 
for VF+dijets at leading order. Also shown for comparison is the equivalent distribution for 
pure dijets. While the VF+dijets distribution is slightly more peaked at zero, there is very little 
difference in the shape of the distributions above a rapidity span of around 1.5, indicating that a 
systematic resummation of radiation between jets is just as relevant for VK+dijets as for dijets. 

The analyses presented above demonstrate that at the LHC, the sub-processes allowing for 
colour octet exchange in the ^-channel dominate the rate for VK+dijets. This observation will 
be the guiding principle for the resummation for PF+jets implemented in HEJ. 



3 W Plus Jets in High Energy Jets 

The framework of High Energy Jets (HEJ) [II 12 constructs explicit approximations to the 
real and virtual corrections to the perturbative hard scattering matrix element at any order. 
Furthermore, matching to the full tree-level results for the first few higher order corrections 
are implemented through a jet merging algorithm. The framework, and the application to the 
description of VK+jets is described in this section. 



3.1 All-Order Amplitudes 

The starting point for the HEJ approach is an approximation to the hard-scattering matrix 
element for 2 — > n partons plus the leptonic decay products of the W . This is built from the 
dominant terms in the High Energy (or Multi-Regge Kinematic) limit, which is defined as: 

Vi,j Sij ->■ oo \p i± \ ~ \p j± \ ^ 
or equivalently i G {1, n — 1} < < y i+ i \pi±\ ~ \pj±\, 

where i,j label the n final state quarks and gluons, ordered in rapidity. There is no constraint 
on the momenta of the e, v. In this limit, the scattering amplitude is dominated by the poles in 
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A y (Forward jet - Backward jet) 



Figure 1: The normalised distribution of the rapidity difference between the two jets, dcr/d(y/ — 
Ub), for W plus two jets at tree-level (black), within the cuts of Eq. ([!]). Also shown for com- 
parison is the same distribution for dijet production, with equivalent cuts on the jets. 




Pi+i 




Pi+l Of V<(q„q 1+1 ) 




Figure 2: Sample kinematics for qQ — > q'g . . .gQ. The diagrams on the right-hand side show 
the dominant kinematic configurations in the High Energy limit where the outgoing partons 
are drawn and numbered in order of increasing rapidity. The effective vertices, V(qi, qi+i), are 



defined in Eq. (10). 
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the i-channel momenta, as depicted in Fig. [2] and defined as 



Qi =PA - (Pe ~Pu) • S — pi - . . . —Pi-l, (3) 

where S = 1 if the W is interpreted as emitted from the quark line with momenta pa,Pi, and 
5 = otherwise. The numbering of momenta is according to decreasing rapidity, A indicates 
the forward moving incoming parton, B the backward moving one. The leading contributions 
to the poles arise from the particle and momentum configurations which allow for a colour 
octet exchange between each neighbouring (in rapidity) pair of partons. We will call these 
FKL configurations. For example, in dijet production ug — > ug and ud — > uggd are FKL 
configurations with that rapidity order, while ug —¥ gu and ud — > ussd are not. Each relaxation 
of a colour ordering induces a relative suppression of 1/s in the squared matrix element and 
hence is subdominant in the High Energy limit. The all-order treatment in HE J currently only 
describes FKL configurations. Other contributions are included order by order with standard 
tree- level matrix elements. 

Scattering amplitudes factorise according to rapidity in the High Energy limit. Based on 
the work of Ref. |11| - |13| we choose a form of amplitudes for each helicity configuration, which 
is a) gauge invariant, b) exact in the MRK limit, and c) sufficiently fast to evaluate that a 
numerical integration over the phase space of the sum over any number of emissions can be 
accurately evaluated. For the case at hand, the tree-level approximation to the square of the 
matrix element is 



.2 /.- 1 \ I ,;i 



2 Tj ■ K >- i^l) (4) 

n {fCA (^ v *(« h * H )'w«'*M>)) 



The bold-face pi indicates the transverse components, with > 0. Each part of this equation 
will now 
currents 



1 1 2 

will now be explained in full detail. \\Sf a f b ->/!/„ Pe P J| is the helicity sum of the contracted 



I c ||2 _ 9\y 



1 



p 2 w -M^ + i T W M W 



2 y \s~, hb T*r~\ 2 , (5) 

/ y 1 ud^fdu e £d 1 ' v ' 

h b ,h2 



where we have picked the case of the extremal partons to represent the scattering ud — > dv^td 
as an example, and the incoming u-quark as the forward-moving parton (i.e. with momentum 
p a ). The spinor contraction for each helicity is simply 

s Jxtd»T~ = jw ^> £ > 1 o) 9 " u uh2 (p^ uhb w ( 6 ) 

with 

V \PW~rPout) (Pin PW ) J 



This is simply the spinor representation of a fermion line emitting a (decaying) W and a (possibly 
off-shell) gluon, as represented by the figure in the equation above. 

In (/^-initiated processes, the W is obviously attached to the quark line. However, for qQ- 
initiated processes, the W emission can be assigned to one or both lines. However, the formalism 
developed needs the VK-emission to be assigned to just one line. Therefore, we assign the W to 
either quark line a or b on an event-by-event basis according to the ratio: 

a: b = e -\Va-vw\ . e -\Vb-vw\ t ( 8 ) 

so that the quark closest in rapidity to the W is favoured, but not chosen exclusively. The event 
weight is adjusted so the outcome is independent of this choice. Picking one specific quark line 
allows for the definite assignment of t-channel momenta in Eq. ^ and the calculation of the 
squared amplitude in Eq. 

For a quark line, the colour factor Kj. is just Cf- We will use the results of Ref. 12 to 
include the case of qg — > q'g ■ ■ ■ gglt by simply using a colour factor depending on the light-cone 
momenta fraction z = Pn/Pb °f the- scattered gluon: 



Kg 



11 + 22 f c A -±-\ + 1 



2 z V " C A C a 



(9) 



It can be seen that K g — > Ca for p n /p b — > 1, in agreement with the well-known result in this 
limit [25 . 

The i-channel propagator denominators are given by ti = qf, with qi defined in Eq. d3J). 
Additional (gluon) emissions are approximated with a string of effective emission vertices f 1 1] , 
given by: 

VP(q i ,q i+1 )=-(q i + q i+l y 

, PA ( <li , Pi+l'PB . Pi+1-Pn\ . 

+ o + + j+PA^Pl nm 

2 \Pi+l'PA PA-PB PA-Pn J (!U) 

P B ( <fi+l , Pi+1 ■ PA , Pi+1 -Pl\ ^ 

- -7T + + -PB^Pn 

2 \Pi+i -p B PB' PA PB-Pl J 

where p g = qi+i — qi is the momentum of the emitted gluon. This simple structure makes it 



very fast to integrate in an efficient phase-space generator [26,27 . This allows the number of 
final state particles to be treated as a variable in the integration. The four-momenta of all final 
state particles is available in every event, allowing arbitrary cuts and analyses. 

An infrared pole will be generated from the phase space integration over the region \p g ^\ — > 
for each i. In these regions, 

-1 4 

— V^isii, qi+l)V^(qi, q i+ i) -»• (11) 

UU+i Pi 

The simple 4/p?-term is therefore used as a real-emission subtraction term (with the integral 
form added to the virtual corrections). Infrared poles are also generated by the inclusion of 
virtual corrections to the t-channel gluon propagator factors. These are accounted for by the 
Lipatov Ansatz [24], which replaces the standard gluon propagator factor 1/ij with 

- -> - exp [a{qi)(yi-i - y { )] (12) 



1 This approach does not include the (kinematically suppressed) effects of quantum interference from the 
emission of W in the cases where there are two incoming quark lines of same flavour) 
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with 

«fe) = VQ^=|'(q 2 /M 2 ) £ - (13) 

We organise the cancellation of the poles from these two sources using dimensional regularisation, 
the subtraction term confined to a small region of real emission phase space with p? < A 2 . 



Further details in Ref. 13 



The final expression for the all-order regularised square of the scattering matrix element is 
then given as 

1 



|-^HEj(fe>Pe,£V})| — ^ J-jyg _ j-j l|5/a/i,->/l/n Pe 



* ** i) ■ [f *>• £ 

n-2 



II [fC A (^^(gi,« + i)^(ffi,«i+i) - p 9 (p 2 < A 2 ))) 



i=i 

n-1 



Y[ exp [uj (qj, \){yj-i - Uj)} , 



0/ \n a s N c 
w fe,A) = log 



7T - A 2 ' 

(14) 



3.2 Merging with Fixed Order 



The regularised matrix elements in Eq. (14) allow for an integration over the momenta of all 



particles. The cross section is then obtained as 



resum _ \^ ST TT { [ Pi± ~°° d 2 Pi± f dj/A f d 3 Pe f d 3 Pu 

W+2J " h k ii \L=o (2-) 3 J 2 ) J ( 27 r)3 2E e J (2^ 2E V 

-2 



h,h n = 2 »=1 



\M^({ Pi , Pe , P „})\~ (15 ) 

p • XajA,fi\%aMa) m %2JBJ 2 [ x b, Qb) 



(2vr) 4 5 2 l^2pk± + Pe± + Pu± 2j ({ Pi }), 



\k=l 



where the first sum is over the flavours fa, fa of incoming partons. In order to stay within the 
relevant kinematics for the formalism, we require the two extremal partons to be part of the 
extremal hard jets identified by the jet clustering algorithm. Furthermore, we require at least two 
such hard jets to be present. These requirements are implemented by the jet observable 2 j in 



Eq. (15). The distribution of any observable can be obtained by simply binning the cross section 



in Eq. (15) in the appropriate variable formed from the explicit momenta. Obviously, multi-jet 
rates can also be calculated by multiplying by further multi-jet observables O^j, O^j, .... 

The simple structure of the HEJ framework makes it feasible to merge with other theoretical 
descriptions, where relevant. By default, HEJ contains matching to fixed-order tree-level ampli- 
tudes in two different ways. Firstly, for flavour and momentum configurations arising in the re- 
summation described above, the approximation to the n-jet production can be reweighted to full 
tree-level accuracy. This is achieved by a merging procedure which clusters all the m momenta 
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generated by the resummation into n on-shell momenta |p^ w ({k})| close to the reconstructed 
jet momenta [113] . The event weight is then adjusted with the ratio of the full n-jet matrix 
element (evaluated by MadGraph |28|) to the approximate one obtained as the a"-expansion of 
Eq. (14) to the n-jet production rate (obtained from Eq. (|4|) 



Currently, reweighting up to 
n = 4-jets is applied. The reweighted resummed cross section is then found as 



oo n 



resum, merged 
T W+2j 



V V TT f r ±= °° ^± I" dyi 

\M^j{{pi,Pe,Pu})\ 2 



d 3 p e 



d?Pv 



s 2 



(2tt) 3 2E e J (2vr) 3 2E V 
XafA,h{ x a,Qa) ■ X 2 fBJ 2 (xb,Qb) 



^ m— jet 

(2tt) 4 5 2 (^2 Pk± + Pe± + PuA ] 2j ({ Pi }), 



(16) 



Kk=l 



with O^j the exclusive m-jet observable and 



w 



n— jet 



M^{[p^{{ Pi }))) 



2 ' 



(17) 



The event sample arising from the equation above is finally supplemented with the kinematic 
configurations (up ton = 4-jets) not arising in this resummation. This is obviously a very naive 



matching, and improvements along a CKKW-L |29 30 procedure should be pursued. Alterna- 



tively, the HEJ resummation could be expanded to cover configurations formally subleading in 
the MRK limit. 

In the meantime, we can assess the importance of the fixed-order contributions with the 
present implementation. Figure [3] shows the rapidity span distribution including different levels 
of matching. Clearly at small rapidity spans, the non-FKL 2-jet matching is dominating the 
cross-section, demonstrating the importance of including these contributions. The impact di- 
minishes as the rapidity span increases (as expected). The further addition of non-FKL 3- and 
4-jet matching (blue and black lines) is visible but less significant for this variable. 

Another important direction for the development of HEJ is work on merging with parton 
showers. HEJ is constructed to approximate the emission of QCD particles at wide angles. 
While this includes emissions with transverse momentum as low as 1 GeV, it does not include 
a description of collinear emissions. This can be added by consistently merging with a parton 
shower program. However, this must be done with some care to avoid double-counting soft 
emissions. A subtraction scheme to merge with the Ariadne parton shower [30] has been 
developed (2l], and further work is ongoing to merge with other parton shower implementations. 

The resummation and merging procedures discussed so far do not rely on a particular choice 
for the factorisation and renormalisation scales. We have implemented four different scale choices 
for fiR = (ip (although they in principle do not need to be set equal) as follows. We list them 
here with the numbering scheme applied in the input file for the program: 
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Figure 3: The rapidity span distribution from HE J for inclusive W plus dijet production for 
the FKL configurations (turquoise), with 2-jet non-FKL added (red), then with 3-jet non-FKL 
matching added (blue) and lastly also including 4-jet non-FKL merging. 

0: Fixed scale of your choice, set in "scale". 

1: Ht/2, where Ht is the transverse sum of all final state particles including p±£± and 
2: P±max ~ the maximum p± of any single jet in each event. 

l/n 



3: The geometric mean of the identified hard jets, ( Yl^ = iP±j 

In addition, there is the option to add logarithmic corrections, mimicking the part of the 
NLL BFKL corrections which are proportional to the LL kernel. These can be included by 
setting the "logcorrect" parameter in the input file to 1. These corrections modify ojo(q,X) of 



Eq. ( 14 ) to give instead 



a s (fj?) C A ( q, \ / a s {fi 2 ) j3 /i 4 




M<1 ^ X)= — ln \ * ) \ 1 + — ln q|A^ J ' (18) 

11 2 

Po =y Nc ~ 3 n f 

while for the real emission vertices the coupling is multiplied by 

a s (/, 2 )^lnp^ 2 ). (19) 

See Ref. [I3j for a full discussion. For the results presented in the next two sections we have 
used fiR = (jLp = Ht/2 and included the logarithmic corrections associated with the scale. This 
last choice is motivated by reducing the impact of the NLL corrections. These choices are not 
an indication of an optimised fit to data, and other scale choices could be studied. 



2 Note that this is strictly greater than the value that can be measured in the detector, where instead the 
momenta of the jets is added (not the original partons). 
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4 Comparison to LHC data 



Having outlined the description of W plus jets in HEJ, in this section we compare the resulting 
predictions to the data we have from the LHC. A recent ATLAS study of W production in 
association with jets pj gave results for many interesting distributions, using the full 2010 data 
sample of 36 pb" 1 . The cuts used in that study, and therefore here, match those of section [2] 
(eq. 0), with the addition of an isolation cut, AR(£,j) > 0.5, applied to all jets. Throughout, 
the HEJ predictions are shown together with a band indicating the variation found when varying 
the renormalisation and factorisation scale by a factor of two in either direction. Obviously, 
this variation is only a rough indication of the true uncertainty of the prediction, with similar 
caveats as for a fixed order calculation. The total time required to generate the event sample of 
this section, and the next, is roughly one day on a single PC. 

We begin in figure [4] by showing comparisons for the distributions of Ht and the transverse 
momentum of the hardest jet, in events of 2,3 and 4 jets. The left plot shows the Ht distribution 
for inclusive W + 2, 3 and 4-jet samples. In all but the bin of lowest Ht for the inclusive Wjj 
production (where the experimental uncertainty is relatively large), the predictions from HEJ 
overlaps with the data, within the quoted uncertainties. It is clear from this plot that at for 
Ht larger than roughly 400 GeV, the suppression from requiring one additional hard jet is only 
roughly a factor 2, and not as- HEJ is developed particularly to deal with the case of a large 
impact from high jet multiplicity. A comparison between HEJ and other theoretical descriptions 
for this distribution has recently appeared in 31 . The right-hand plot of figure [4] shows the 



transverse momentum distribution of the hardest jet in each event, again for inclusive W + 2,3 
and 4-jet events. This, more inclusive, variable is less sensitive to the description of additional 
radiation like that included in HEJ, and still one sees good agreement between the HEJ prediction 
and data. 

Figure [5] shows the distributions for the invariant mass of the hardest 2, 3 and 4-jets in the 
event. This is more sensitive to the topology of the events, in addition to the overall momentum 
scale. As the number of jets increases, the peak of this distribution moves away from the 
kinematic minimum to higher values of invariant mass. Over this wider range of momentum, 
we again find a very close agreement between the HEJ predictions and the data. 

The final plots we show in this section, in figure [6j probe the relative position of the two 
hardest jets in each event. Firstly, the left-hand plot shows the distribution of the difference 
in azimuthal angle between the two hardest jets. This is peaked at A<f> = tt. For pure dijet 
production at tree-level, the azimuthal distribution is a delta-functional at tt, and higher order 
corrections smear out this distribution, which however remains peaked at Acf) = tt. The fact 
that data for Wjj shows a similar structure for the azimuthal distribution indicates that this 
process proceeds less like one parton recoiling against a W, and then splitting into further jets 
and more like a dijet scattering with a VF-emission. This indeed is the mechanism implemented 
in HEJ, and it describes the data well across the distribution. 

The right-hand plot in figure [6] shows the distribution of the difference in rapidity between 
the two hardest jets. It is clear that the HEJ description is slightly high at the peak around 
Ay = 0. This is precisely the region where the impact of matching to fixed-order is largest, and 
the impact of the naive scheme currently applied will be greatest. This is an important area of 
future development for HEJ. However, we see that across the analyses, the description of VF+jets 
as currently implemented in HEJ describes existing LHC data very well. 
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Figure 4: Left: The W + jets cross section as a function of Ht, the scalar sum of the transverse 
momenta of the jets, the charged lepton and the missing transverse momentum. Right: The 
W + jets cross section as a function of the transverse momentum of the hardest jet in the event. 
The data points, in this and subsequent plots in this section, are taken from Ref. [2j. 
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Figure 5: W + jets cross section as a function of the invariant mass of two, three and four jets 
separately. 




Figure 6: Left: The W + jets cross section shown as a function of the absolute difference in 
azimuthal angle between the two hardest jets, \4>j e ti — 4>jet2\- Right: The W + jets cross section 
as a function of rapidity difference of the two hardest jets. 
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5 Probing Higher Order Corrections 



We have seen that the HEJ approach gives a good description of the current LHC data for 
the production of a W in association with jets. In this section, we now turn our attention to 
variables which may better distinguish between the predictions obtained in standard approaches 
(like NLO, a parton shower, or a combination thereof), and that implemented in HEJ. 

Throughout this section, the HEJ prediction is shown as the black central line. The result 
of varying the renormalisation and factorisation scale by a factor of two in either direction is 
shown by the cyan band. The central scale choice chosen in this section is p±ma.x, the maximum 
transverse momentum of any of the jets. 

Since all of the remaining studies are of ratios of cross sections with a correlation between 
numerator and denominator, the statistical uncertainty cannot be estimated by the usual prop- 
agation of errors derived by assuming no correlation. Instead, the Monte Carlo uncertainty was 
evaluated by splitting our generated events into 12 samples of equal size, and calculating the 
distributions from each of the 220 possible ways of selecting 9 out of these 12 samples. The 
statistical error band is then defined such that 68% of these predictions lie within the central 
line obtained from all 12 samples. 

The first observable we will consider is the average number of observed jets. The HEJ 
predictions for this are shown in figure [7] as a function of Ay (left) and Ht (right). The top 
plots show results for inclusive Il^+dijet samples, but the predicted average number of jets in 
the events rises to around 3 in each case, emphasising the importance of terms at higher orders 
in a s for a reliable jet count. This may play an important role in discriminating SM vs. BSM 
contributions to the same channel. Both of the regions studied here are important regions of 
phase space. This variable has been studied by the ATLAS collaboration in the context of dijet 
production |14| , where significant effects beyond fixed order were seen. 

In the bottom row of figure [7J the same variable is plotted but now restricted to events 
with three resolved jets or more. Again, we see that higher orders contribute significantly here, 
with the average number of hard jets reaching around 3.3 and 3.5, when shown as a function 
of Ay (the rapidity difference between the most forward and most backward hard jet) and Ht 
respectively. 

In Ref. [3l], predictions for the average number of jets in W+dijet events were compared 
between four theoretical approaches: HEJ, a pure Wjj NLO calculation [4j[6j[7], the "NLO exclu- 
sive sums" approach (a brute force method to combine NLO calculations of different orders) and 
the Sherpa 32,33 MEPS (34-36 scheme which combines tree- level matrix elements of different 



orders in a s with a truncated parton shower. Large differences were seen in the predictions 
between these theoretical descriptions and an experimental study of this variable would further 
our understanding of the nature of QCD radiation in the high-energy environment of the LHC. 

The ratios of the inclusive jet rates are of course slightly less sensitive to additional radiation 
than the average number of jets. In Fig. [8] (top row) we plot the 3-jet to 2-jet rates, along with 
the ratio of the inclusive 4-jet to 3-jet rate (bottom row). Once again it is clear that the impact 
from higher orders in large. For illustration in the top row, the ratio between the tree-level 
3-jet and 2-jet rates has also been plotted. Although this contains no systematic resummation 
of higher orders, in fact the leading-order result rises higher than that of HEJ. The ratio of the 
inclusive 3-jet to 2-jet rate was also studied in Ref. 31 for the different theoretical descriptions 
listed above. As expected, smaller differences were seen between the approaches than in the 
predictions of the average number of jets, but even so, the experimental data could probably 
select a preferred description of the inclusive jet rates, especially at large Ht- 

The ratio of the inclusive 4-jet to 3-jet rate is also shown in the bottom row of figure [8l The 
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Figure 7: Top Left: The average number of jets as a function of the rapidity difference between 
the most forward and most backward jet and Top Right: the average number of jets as a function 
of Ht- Bottom Row: as top row, except restricted to events with 3 jets or more. 
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Figure 8: Top Left: The ratio between the inclusive 3-jet and 2-jet rates as a function of the 
difference in rapidity between the most forward and most backward jet; Top Right: The ratio of 
the inclusive 3-jet to 2-jet rate as function of Ht- The dotted line in each plot shows the ratio 
between the 3j and 2j tree-level calculations as a function of Ay and Ht respectively. Bottom 
Row: as top row, except the ratios are now of the inclusive 4-jet rate divided by the inclusive 
3-jet rate. 
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increase here is somewhat smaller than that seen for the ratio between the 3-jet and 2-jet-rates, 
but it still rises to 30% and 45% as a function of Ay and Ht respectively. 

Finally, in Fig. [9] we compare the predictions from HE J for the average number of jets vs. Ay 
for inclusive VF+dijets at the 7 TeV and the 8 TeV LHC. The result for 8 TeV shows only a 
very modest increase in the average number of jets. 

In this section the predictions from HE J for the average number of jets and ratios of inclusive 
jet rates have been shown. These show a large degree of sensitivity to additional hard QCD 
emissions beyond the leading order. These higher-order effects will only increase with the centre- 
of-mass energy of the LHC collisions. 

6 Conclusions 

We have described the application of the High Energy Jets (HEJ) framework to the production 
of a W boson in association with at least two jets. HEJ resums systematically the contribution 
from multiple hard emissions (including also the leading virtual corrections). The process of 
W+dijets offers a key testing ground for our understanding of the behaviour of the Standard 
Model at the LHC, and will in turn be important for many directions including analyses of Higgs 
boson couplings and searches for new physics. 

The predictions of HEJ were seen to give a very good description of the distributions studied 
with the 2010 data set. We further considered observables which are designed to be sensitive to 
the final state configuration of hard jets, and thus probe the perturbative description. We saw 
that the impact of higher orders on the inclusive dijet sample are large, for both the average 
number of jets and the ratios of inclusive jet rates in various regions of phase space. This offers 
the possibility of observing directly the impact of the BFKL-inspired resummation offered by 
HEJ. 

The implementation of the formalism in the form of a fully flexible partonic Monte Carlo, 
can be downloaded at http://cern.ch/hej. 

Higher order QCD effects have already been observed in data for pure jet production, and 
we look forward to future LHC analyses of W+jet production to further our understanding of 
physics at these new energy scales. 
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